Decoherence, Chaos, and the Correspondence Principle 
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We present evidence that decoherence can produce a smooth quantum-to-classical transition 
in nonlinear dynamical systems. High-resolution tracking of quantum and classical evolutions 
reveals differences in expectation values of corresponding observables. Solutions of master 
equations demonstrate that decoherence destroys quantum interference in Wigner distribu- 
tions and washes out fine structure in classical distributions bringing the two closer together. 
Correspondence between quantum and classical expectation values is also re-established. 
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The status of quantum-classical correspondence for dy- 
namical systems is somewhat unclear and perhaps even 
controversial. In a nonlinear system, a single classical tra- 
jectory initially centered on a quantum wavepacket can 
quickly diverge from the motion of the centroid given by 
the quantum expectation values of x and p (this defines 
the Ehrenfest time, c.f. Ballentine et aim Ref. Q). How- 
ever, when the single trajectory is replaced by a classical 
probability distribution to give the correspondence prin- 
ciple a better chance, the situation is far less certain. 
Thus, some authors have reported evidence for a break- 
down of the quantum-classical correspondence in chaotic 
systems || , while others have argued that it can be pre- 
served when stated in terms of the expectation values of 
dynamical variables Q. In yet another line of investi- 
gation it has been shown that, even in chaotic systems, 
semiclassical methods are successful for times longer than 
was previously believed possible || . 

In integrable systems, a rapid divergence between clas- 
sical and quantum evolutions can occur for initial condi- 
tions near an unstable point: This can force the system to 
undergo a "double-slit experiment," yielding very differ- 
ent outcomes in the two cases. This breakdown of corre- 
spondence is caused by the coherent interference of frag- 
ments of the wavepacket and occurs on a short timescale, 
on the order of the system dynamical time. However, 
for generic initial conditions, quantum and classical ex- 
pectation values are expected to diverge on a timescale 
inversely proportional to some power of h |@J|. This is 
regarded as sufficiently slow to cause no difficulties with 
the classical limit of quantum theory. 

The first purpose of this Letter is to investigate the 
quantum-classical correspondence at the level of expecta- 
tion values for chaotic dynamical systems, where a break- 
down may be anticipated on a much smaller logarithmic 
timescale t% ~ ln(C/7i) H,|||3)- We show below that such 
a loss of correspondence does occur though its magni- 
tude is not in conflict with earlier results Q|. While the 
timescale on which the violation occurs varies with the 



particular initial condition chosen, it is consistent with 
the logarithmic time ta- 
in order to affect a quantum-classical comparison we 
chose Gaussian packets as initial states, with positive def- 
inite Wigner distribution functions. The time evolution 
was then performed using either classical or quantum dy- 
namical equations. In the chaotic systems investigated 
here, for generic initial conditions (i.e., Gaussian packets 
randomly sampling the chaotic part of the phase space) , 
differences between quantum and classical expectation 
values stay small for some time and then abruptly in- 
crease. After this divergence time, the differences remain 
modest, typically ~ 5 - 10% (for h = 0.1). The Wigner 
function begins to differ considerably from the classical 
phase space distribution at relatively early times 0] . 

Our second purpose is to show that the discrepancy 
between quantum and classical evolutions is drastically 
decreased by even a small coupling to the environment, 
which in the quantum case leads to decoherence ||. 
There are two limiting situations in the study of decoher- 
ence in dynamical systems. In the first case, special ini- 
tial conditions (such as Schrodinger cat states) are used 
to study the destruction of interference already present in 
the initial state but with simple system dynamics, typ- 
ically taken to be linear However, since quantum 
interference is dynamically generated only in nonlinear 
systems |]lo| , the competition between generation and de- 
struction of quantum coherence cannot be investigated. 
The results reported in this Letter are from a study of the 
second type, where the system dynamics is nonlinear, but 
the choice of initial states is kept deliberately simple so 
as to focus only on the role of dynamically induced inter- 
ference (as distinct from that present in the initial state). 
We show below that differences at the level of expecta- 
tion values are sharply reduced due to decoherence and 
the effect on correspondence in phase space is even more 
spectacular. In the parameter regime investigated here - 
essentially the border between quantum and classical - 
some remnants of quantum coherence may still survive. 
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However, such small-scale coherence has apparently little 
effect on the correspondence of the expectation values. 

A mechanism responsible for the quantum to classical 
transition should explain not just how expectation values 
can converge to the same answer, but also lead to com- 
patible effective phase space distributions. A common 
approach is an appeal to coarse-graining, a formal proce- 
dure implemented typically by convolving the individual 
distributions with a Gaussian and then comparing the 
two resulting coarse-grained distributions Q. This ap- 
proach has three defects. First, as a formal mathematical 
procedure it can always be inverted, and thus offers no 
physical insight. Second, this coarse-graining does not 
alter the dynamics, and hence cannot improve the con- 
vergence of expectation values. Third, for the classical 
system, the notion of a trajectory is lost and along with 
it the notion of Lyapunov exponent. 

In contrast to the coarse-graining approach, decoher- 
ence provides a dynamical explanation of the quantum 
to classical transition by taking into account interactions 
with an (external or internal) environment of the sys- 
tem - degrees of freedom that effectively monitor and, 
therefore, select certain stable or "pointer" observables 
destined to become the classical variables [[|[ll]. The 
simplest models of this type lead to master equations for 
the reduced density matrix for the system [Q. Diffu- 
sion terms in these equations automatically coarse-grain 
the distributions, now a physical effect of the coupling to 
the environment rather than a mathematical trick. The 
degree of coarse-graining is determined by the interplay 
between the dynamics of the system and the nature and 
strength of the coupling with the environment. More- 
over, the effectively classical master equations that de- 
scribe the post-decoherence dynamics admit a Langevin 
description of trajectories allowing for the existence of a 
Lyapunov exponent. We demonstrate below that deco- 
herence dramatically improves the correspondence of the 
expectation values, leads to the existence of a single effec- 
tive phase space distribution, and allows for a Lyapunov 
exponent to exist at late times. All the deficiencies of the 
coarse-graining approach are therefore overcome. 

We restrict attention to bounded, one dimensional, 
driven systems. Tools employed are very high-resolution 
simulations of the time-dependent Schrodinger, quantum 
and classical Liouville, and master equations recently im- 
plemented on massively parallel computers Jl| . The nu- 
merical results discussed below are for the driven system 
considered in Ref. H], with Hamiltonian, 

H = p 2 /2m + Bx 4 - Ax 2 + Axcos(ujt) . (1) 

We used a parameter regime (m = 1, B = 0.5, A = 10, 
A = 10, uj = 6.07) in which a substantial area of phase 
space is predominantly stochastic, with the finite-time 
Lyapunov exponent A ~ 0.4 — 0.5. Gaussian phase space 
distributions, typically minimum uncertainty wavepack- 
ets, were employed as initial conditions to sample the 



evolution in the stochastic region. Other Hamiltonians 
studied include the driven Duffing system and a two- 
dimensional, truncated Toda potential. These systems 
yielded similar results; a detailed presentation will be 
given elsewhere Q. 

The specific model of decoherence used here is the weak 
coupling, high temperature limit of quantum Brownian 
motion Jl^ . In this limit, dissipation can be ignored 
at early times, and only the diffusive contributions in 
the master equations |||| need be kept. The diffusion 
constant was chosen to be small enough so that, over the 
timescales of interest, the energy increase was negligible, 
and changes in the classical phase space structure were 
only perturbative. For the Hamiltonian ([l]), the quantum 
master equation in terms of the Wigner function is 
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where dV/dx = 4Bx 3 — 2 Ax + Acos(wi) and L q is 
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The classical (L q = 0) Fokkcr-Planck equation is 
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This dynamics is equivalent to the Langevin equation 
mx = —V'(x) + F(t) where F(t) denotes a Gaussian, 
white noise. The quantum Schrodinger and master equa- 
tions, as well as the classical Fokker-Planck equations 
were solved using a high-resolution spectral algorithm 
pUf . In the absence of diffusion, the classical Liouville 
equation was solved as an N-body problem, the distribu- 
tion being sampled by at least ~ I0 5 particles. Numeri- 
cal checks included carrying out simulations at different 
spatial and temporal resolutions and a direct verification 
for the moments obtained from the codes by substitut- 
ing them in the BBGKY-like moment evolution hierar- 
chy and verifying that this set of equations is satisfied to 
a relatively high order. Lyapunov exponents were com- 
puted using the techniques described in Refs. Jig ]. It 
was verified that at the noise levels used, the finite-time 
Lyapunov exponents from the Langevin equation agreed 
with those computed from the Hamiltonian dynamics. 

The exponential instability characteristic of chaos 
forces the system to rapidly explore large areas of phase 
space and to interfere on a timescale set by when the 
wavefunction has spread over much of the available space 
and the Moyal corrections arising from its nonlo- 
cality have become comparable to the classical force H: 



t n ~ X-'HxSp/h) 



(4) 



where A is the Lyapunov exponent, Sp is the mea- 
sure of dispersion in the initial conditions, and x — 
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\Z\(d x V/d xxx V)\ is a measure of the nonlinearity in the 
potential averaged over the accessible space. We inves- 
tigated the evolution of several expectation values, such 
as (x), higher order moments such as ((x — (x)) n ) (with 
a maximum n = 4) , and expectation values of variables 
that in principle include all moments. Some of our re- 
sults are shown in Fig. 1. In all cases, and for all of 
the investigated initial conditions we found good agree- 
ment between the quantum and classical results during 
the initial portion of the evolution. (This initial period 
was longer than the Ehrenfest time, consistent with the 
results of Ballentine et al The onset of the discrep- 
ancy depended on the initial condition, but in all cases 
was a factor of a few larger than the dynamical time. The 
discrepancy saturates to typically no more than 10% of 
the expectation values thereafter. The value of h was var- 
ied to test for logarithmic scaling: while the results are 
consistent with (0), the dynamic range of the simulations 
is insufficient to make a more precise statement. 
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FIG. 1. Classical and quantum expectation values (as) 
as a function of time. The initial condition is a minimum 
uncertainty Gaussian wave packet with (a;) = —3, (p) = 
8, ((as - {x}) 2 } = 0.0025, and ((p - (p)) 2 ) = 1. Ti is set to 
0.1. This yields tn ~ 4. The central vertical bar denotes the 
average divergence time and the left and right bars, the min- 
imum and maximum respectively, for ten initial conditions 
that randomly sampled the chaotic phase space. 

We have therefore good evidence that in isolated 
chaotic systems, the quantum-classical correspondence 
defined at the level of expectation values is lost rela- 
tively quickly due to dynamically generated quantum in- 
terference. This is best appreciated by comparing classi- 
cal phase space densities with quantum Wigner distribu- 
tions. As the wavepacket spreads and folds, the Wigner 
function becomes dominated by small scale interference 
which saturates on a scale set by the size of the system 
in both momentum and position: Sp = h/L, 5x = h/P. 
Eventually, the Wigner function is unable to track even 



the "backbone" of the classical phase space distribution 
and becomes a complicated looking interference pattern 
in which the classical phase space structure can no longer 
be distinguished [Fig. 2(a)]. Fine scale structure in the 
interference pattern (i.e., oscillations within an h box) is 
clearly apparent. On the timescales probed in our simu- 
lations, the Wigner function has not reached the stage of 
"structure saturation," i.e., smoothness on a scale ~ h 
in phase space || . 

Our results are encapsulated in Figs. 2. As can be 
seen from comparing the decohered Wigner function [Fig. 
2(b)] with the classical distribution given by the solution 
to Eq. (H) [Fig. 2(c)], decoherence markedly improves 
the correspondence at the level of distribution functions, 
radically changing the unitarily evolved Wigner function 
of Fig. 2(a) by effectively smoothing it over scales ||; 
Ap ~ o~ c = ^/2Z?/A in momentum (which is translated 
by dynamics into a smoothing in position). In our case 
o~ c — 0.3 and \ ~ 0-6, which implies that we are on 
the border between quantum and classical regimes (de- 
fined respectively by whether a c x is large or small com- 
pared to h H). Even though the distributions in Figs. 
2(b) and 2(c) are very close, the Wigner function still 
contains traces of local quantum interference. However, 
this makes only a minor difference to expectation values, 
and tends to vanish as the evolution proceeds. We note 
that with our choice of parameters, the diffusion term 
affects the evolution of classical and quantum expecta- 
tion values roughly to the same extent (Fig. 1). Deco- 
herence destroys the interference pattern in the Wigner 
function, while at the same time, noise smooths out the 
fine structure of the classical distribution in such a way 
that quantum and classical distributions and expectation 
values both converge to each other. Thus, one concludes 
that the decohered quantum evolution does go over to 
the classical Fokker-Planck limit. 

In summary, we have provided evidence that in a quan- 
tized classically chaotic system, for fixed ft, classical and 
quantum expectation values diverge from each other af- 
ter a time ~ t%. In the case studied here the discrepancy 
is < 10% of the typical magnitude of the expectation val- 
ues. Decoherence was shown to substantially reduce this 
discrepancy as well as to bring the Wigner and classi- 
cal distributions very close to each other. In this com- 
bined sense, decoherence restores the quantum-classical 
correspondence. Our results complement previous stud- 
ies which have focused more on the phase space aspects 
of the correspondence and the destruction of dynamical 
localization by noise and dissipation [jfffl. 
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FIG. 2. (a) Wigner distribution function from a solution 
of Eq. ^ at time t = 8T, where T is the period of the driving 
force. The diffusion constant D — 0. The box represents a 
phase space area of Ah. (b) Wigner distribution function at 
time t — 8T, with diffusion constant D = 0.025, illustrating 
the destruction of large scale quantum coherence, (c) Classi- 
cal distribution function from a solution of Eq. Q at time 
t = 8T, with diffusion constant D = 0.025. 
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